% MT fwd & inverse computaton - keller & kaufman 
% latest date 04.01.2002
% inverse problems ----(???)
% fwd nlayer bug fixed 02.1.2002

% t = period (input - HZ)
% rho = layer resisitivity ohm-m
% h = layer thickness (IN KM)
% j9XY;

mu = 4*pi*10^-7;
d1 = sqrt(1/mu);
%h   = [d1];
%rho = [1,0.01];
h = [200];
rho = [20,2000];

frq = 10.^[-4:.1:4];

t = 1./frq;
[rh,ph]=mtfw(rho,h,t);
subplot(211)
loglog(t,rh,'r.');
subplot(212)
semilogx(t,ph,'r.');



